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^^ ' Abstract 

<N . 

We consider two dimensional real-valued analytic potentials for the Schrodinger equation which are 
j^ . periodic over a lattice L. Under certain assumptions on the form of the potential and the lattice L, we 

can show there is a large class of analytic potentials which are Floquet rigid and dense in the set of 
C°°(R^/L) potentials. The result extends the work of Eskin et. al, in "On isospectral periodic potentials 
(T^ . in R", 11." 

X . 1 Introduction 

■^ , The focus of this paper is the class of Schrodinger operators 



P : u{x) ^ (-A + q{x))u{x), 
where 



2 q2 



i=i •' 
and 

O ' q{x) : M^ ^ R 

PsJ , is a real- valued periodic potential over a lattice, L C M . In other words we have 

^— V I q{x + d) = q{x) \fd G L. 

(N ■ 

We will study the question of spectral rigidity for the operator P and derive results which could extend to 
R" for n > 3. We consider the set of A in M for which the self-adjoint eigenvalue problem 

^1 Pu{x) = Xu{x) u{x + d) = exp {2Trik ■ d) u{x) (1.1) 

H ■ 

Cu ' has a solution for k in R^ and d in L. When there is a nonzero solution to (|l.ll) we say that A is in 

Specfe(— A -I- q). We refer to 

IJSpecJ-A + g) 

feeK 

as the Floquet spectrum. However, when k = 0, wc simply say 'spectrum' which we denote by Spec(— A + q). 
Two potentials q and q are Floquet isospectral if 

SpeCfe {-A + q) = Specfe {-A + q) Vfc e R 

and isospectral if Spec(— A -I- q) =Spec(— A -I- q). Following the convention in [3], we consider a potential to 
be Floquet (spectrally) rigid if there are only a finite number of potentials modulo translations which are 
Floquet isospectral (resp. isospectral) to it. 

In [2] , Eskin et al. showed that under the assumptions 



1. q is real analytic 

2. L has the property |d| = \d'\ => d = ±d' for all d, d' in L 

then Spcc(— A + q) determines SpeCfc(— A + q) for all k in R. 

In the sequel to [5], [3], Eskin, et al., show that there is a set of analytic potentials satisfying the 
conditions (1) and (2) which are dense in C°°(M^/L) such that if q{x) is in this set, then q{x) is Floquet 
rigid. Furthermore, there is a smaller, but still dense set of analytic potentials in C°°(M^/L) such that if 
qix) is in this set and q{x) is Floquet isospectral to q{x) then, q{x) = ^(ia; + a) where a is an arbitrary 
constant. Under the assumptions (1) and (2), if a potential in M^ is spectrally rigid (resp. unique) then it 
is Floquet rigid (resp unique), so their results arc also true with the words "Floquet rigid" (resp. unique) 
replacing " isospectrally rigid" (resp unique). The main result of this paper is to show that there is a more 
general class of potentials which satisfy the conditions for Floquet rigidity than in [3] . 

2 The Isospectral Manifold in M^ 

In M} the structure of the isospectral sets of periodic potentials has been well studied and contains many 
results which are useful in higher dimensions. In M^ the Schrodinger operator becomes Hill's operator. 

where q{s) has period 1 and is real-valued. We start by assuming that q is at least three times differentiable, 
so that wc can use many of the standard results which may be found in Magnus and Winkler, [B]. For the 
rest of this paper, wc will also assume that q{x) has mean zero. We look at the set of A where there is a 
solution to 

-^^ + q{^ms) = M{s) (2.1) 

</>(s + l) = (-l)™</>(s). 

The scalars A are known as the periodic and anti-periodic eigenvalues. Through curious use of notation, the 
scalar, A^, denotes the eigenvalue corresponding to the eigenfunction (j)^{s + 1) = (— l)™(?!)^(s) so that 

Ao < Aj" < A+ < A2" < A+... (2.2) 

Hence the periodic spectrum consists of {A^, m even} and the antiperiodic spectrum is {A^, m odd}. 
If we change the problem (|2.1I) so that (f){s) obeys the boundary condition 

0(0) = 0(1)^0, 

then the associated spectrum is called the Dirichlet spectrum. The Dirichlct spectrum arc denoted fJ-mil) 
and they interlace the periodic and anti-periodic spectra. We will often use the fact 

|A+-A+|=0(|m2-n2|), (2.3) 

and find it worthwhile to mention it here. Although A^ < )^m+i^ i^ i^ possible to have A~ = A^. The 
spectrum of 

d^ 



as an operator in L^(K) is 



U I- ™' "^m+l-l 



m=0 



Each of the intervals [A+ , A„_|_]^] in the union above is called a "band" , or interval of stability. The complement 
of the set of bands is union of the intervals (A~ , A+ ) which are called "gaps" or intervals of stability. In each 
gap, the operator — ^ +9(5) does not have a bounded eigenfunction. A gap is referred to as open whenever 
A,7j < A+ and closed if A,7j = A+ . The length of a gap is denoted as 7„i . 

In [1] Garnett and Trubowitz gave a compete characterization of the gaps for q in ^^[0, 1]. 

Theorem 1. 01/ Let Jn? n > 1, be any sequence of nonnegative numbers satisfying 

TI>1 

Then there is a way of placing the sequence of open tiles of lengths 7^7 f^ > 1 ^ order on the positive axis 
(0,00) so that the complement is the set of bands for a function q in L^[0, 1]. In other words, the map 

q ^ i{q) -^ {in{q)}n>i, (2.4) 

from L^[0, 1] to (P)^ , is onto. 

Furthermore if we multiply the gap lengths 7^, by e where e is in [0, 1] then the map (|2.4p is still onto. 
The fundamental result in M} is that the set of analytic periodic potentials M(e) with the same periodic 
and anti-periodic spectra is equivalent to a torus with dimension equal to / [7] . Here / is the number of m 
for which A~ < A^. The coordinates am{q), on this manifold with m referring to the rri*'' gap on q{s), are 
related to the Dirichlet spectra and the gap lengths. They are defined as follows 

■ 2 f \ Mm(g) - A„ n TT 

'^m ^ Am ^ ^ 

where fim{q) is the Dirichlet eigenvalue for q such that A^ < fJ^miq) < A^^. These coordinates arc further 
discussed in Section 4. 

Finally we will need the fact that all the gap lengths are exponentially decreasing if and only if q{s) is 
real analytic. Whenever q has only a finite number of open gaps, then q must be real analytic, [5]. The 
analyticity of q{s) with finitely many gaps is crucial in many of the proofs of the theorems in this paper. 

3 Review of Necessary Results in M" 

Wc outline some necessary results and definitions from [2] and [3] which will be used in the rest of this paper. 
Let L be an n-dimensional lattice generated by n vectors i'i,f2, ■■■,Vn- We can then consider it's dual L* 
where 

L* = {(5e]R" -.S-vezyveL}, 

to be generated by some basis 61,62, ■■■,Sn. A function is periodic over the lattice L if q(x + d) = q{x) for 
all d in L. For any arbitrary lattice L satisfying condition (2) and basis fixed as above, let §* be the set of 
fundamental directions for L, that is 

§* = {(5 e L* : (5 • d = 1 for some d eh}. 

It is clear that whenever 6 is in §* then —S is also in this set, so we reduce the set to S by only picking 6 in 
§*. Therefore any element of L*/{0} has a unique representation as mS with 6 in S and m in Z. 
If (7 is a function which is periodic over L, then it has the following Fourier series representation 

q{x) = 2, ^<5 cxp {2ni6 ■ x) 

with 

as = , , / q{x) exp {-2Tri6 ■ x) dx 
r 



where F the fundamental domain of the lattice L as given by 

r = {siWi + ... + s„w„:0<s, < 1}. 



If we write 



then we have that 



|(5pg<s(s) = y^ ans exp(27rws) 



q{x) = 2^ 2_^ '^"'5 '^''^P (27''Jn(5 • x) = 2^ \5\^<ls{5 ■ x) 

where each qs{s) is a periodic potential on R^. These one-dimensional potentials qs{s)' s are called directional 
potentials. The assumption that q{x) has mean zero is equivalent to setting ao = for all the directional 
potentials. 

Theorem 2 in ([2], [3]) states that 



Theorem 2. S'pec(— A + q) determines 






Spec^ [--^ + qs{s)] ySeS,ke 



The theorems in M^ we mentioned will help reduce the study of periodic potentials in M" to the study of 
IS} potentials, about which much more is known. 

4 Potentials in R^ 

Following j3], for the rest of this paper we assume that the elements of the lattice L satisfy condition (2) as 
stated in the introduction, and we consider analytic periodic potentials q{x) such that q{x + d) = q{x) for all 
d in L. We also only consider potentials with a finite number of directional potentials. For this section, we 
make the additional assumptions that the number of gaps in each direction Sj is finite, and that there are 
at least 3 directions. This setup differs from [3] where two of the directional potentials were fixed translates 
of the one gap potentials and the other directions were viewed as perturbations of the zero potential. 
Under these assumptions we can simplify the form of q{x) as follows 

s 
q{x)^Y.\^,\\,{S,-x). (4.1) 

Each one dimensional directional potential qj {Sj ■ x) corresponds to a one dimensional operator with corre- 
sponding eigenvalue and eigenfunction pair (A, (j){s)) satisfying 

-^Hs) + q,i.sMs) = Xcj,{s). (4.2) 

In order to simplify the computations needed in this paper we make the following assumptions (*) 

1. S^ = Si + S2 

2. qi,q2 and 53 have the same number of open gaps 

We will discuss how, given sufficient time and energy, using spectral invariants and the standard perturbation 
techniques that one could remove the assumptions (*). The invariants are derived from the trace theorems. 
If we let the fundamental solution of the heat equation 

du 

— = Au-qu u{0,x)^f{x) (4.3) 



onR" be G{x,y,t) then 

^ exp(-A<) = ^ exp(-27rifc • d) / G{x + d, x, t) dx (4.4) 



Therefore if one knows Specfe(— A + q) for all fc, then one knows 

I G{x + d,x,t)dx Vt>0,deL (4.5) 

r 

In [2] and [3], they derive Theorem [5] from the asymptotics of 

/ G[x + Nd + e, X, t) dx Vi > 0, d e L (4.6) 

r 

as A^ — ^ oo. 

Theorem [5] has the consequence that the set of real-analytic q{x) isospectral to q{x) can be identified 
with a subset of a real analytic manifold 

A/ = Ti X Tz X ... X Ts- 

Here each torus Tj has dimension equal to the number of open gaps associated to each directional potential 
q{6j ■ x); we call this set Ij. This manifold M has dimension X) l-(jl ~ ^- Again, the coordinates on the 

j 
manifold aj^rnil) are given for each j by (|2.5p . 

In our case, we would like our set of potentials which we will call M{e) to have open gap lengths which 
are parametrized as follows. Let Eq denote the set 

{(j,m):(j,m) = (l,l),(2,l)}, 

and El denote the set 

{{j,m) : j < 2,171 > 1}. 

Now we let e be the vector with four components (ei, £27 £3, £4) so we can parametrize the new gap lengths 
so they depend on e and 7 as follows 

7j,m(e,7) = ej7j,m for (j, to) G Eq 

7i,m(e,7) = (-4:1], ni for (j,to) G Ei 

73,m(e,7) = e373,m for m G I3 

lj.in{e-i 7) = e47i,m for j > 3, jti G Ij 

and are associated with the potential g(e, x, a). Here, suppressing the "g" , we have a — {aj,m\ is the rescaled 
vector of coordinates, where for each directional potential, the coordinates arc given by (j2.5p . Notice that 
we have also written our gap lengths in terms of finitely many parameters and this does not destroy the fact 
the mapping (|2.4p is onto and in this case analytic. 

The following spectral invariants are derived from higher order terms in the asymptotics of 14.61 in [3] 
which we will use in our computations: 

Theorem 3. The periodic and anti-periodic spectra for the one dimensional potentials qs{x) which form 
qix) and the invariants 

$5,,m(e,a) (4.7) 

= $j,,„(e,a)= / |/i(e,a;,a)p(0^„(e,(5j •a;,a))^da; 



when ^jm > ^j m '^'^'^ 



*d-,,,»(e,a) = $j-™(e,a) (4.8) 

\h{e,x,a)\^ ((0+„(e,(5j ■ x,a))^ + {(t)-^{e,6j ■x,a))^) dx 



when X\^ ~ ^jm '^laj/foe recovered from the spectra of q{x). Here a ~ {ctj.m] is the collection of coordinates 
associated to each gap length and we have set 

/i(e, a;, a)= ^ — —qe{e,e-x,a) 






^3 



with Sj ■ dj = 0, and dj of minimal length. 

Setting ^g m(^''-^) ~ ^j.mi^iCt), then the number of invariants with A+ > A~ has dimension equal to 
the manifold M{e). We would like to show that the Jacobian determinant of the invariants with respect to 
the coordinates a is nonzero so that we may apply the implicit function theorem. 

We will primarily be calculating the spectral invariants for potentials at a specific parameter e = eq. We 
let eg be the vector with (ei, 62, 0, 0) where ei and 62 are in (0, 1). When e = eq the potential q{eo, x, a) has 

7i,m(eo,7) = ej7j,m for(j,m) £ Eq 
7j,m(eo, 7) = for {j, m) e £:g 

for gap lengths. The potential q{eQ,x,a) is therefore the sum of 2 potentials with only one gap, one in 
each direction Sj, j = 1,2. The rest of the directional potentials are zero. While the limit q{eo,x,a) 
coincides with the form of the potential as calculated in [3], one specific difference remains- the first two 
directional have finitely many gaps, they are not just translates of the p function. We will Taylor expand 
the Jacobian determinant with respect to £3 around e 7^ eo and use these computations to show that the 
Jacobian determinant for certain fixed a is not identically zero. 

For the rest of this paper, we let p(s+^, r) denote a general normalized Weierstrass p function. Whenever 
the parameter r is real and greater than zero, then p{s + ^,r) is real- valued with periods 1 and r [5]. The 
real- valued p- function is always even about i, and by a theorem of Hochstadt [S], all one gap potentials 

are translates of the p-function. The directional potential, in the limit, qj{£o, s, a) = p{s + ^ + vj,tj) has 
eigenfunctions which satisfy the following equation: 

--7^0(eo,s, a) + gj(s)0(eo,s, a) = A^(eo,s,a). 
where qj{eoj 0, a) = p(-j- + Vj, Tj) has bands given by 

[-^^ (1) '-^ (^)] ^ [-^ (t) '+^)- ^^-^^ 

Aligning the classical elliptic function theory with spectral theory [1] we have that, 

1 \ ^ f iTj + l\ ^ _ f ITj 



2) ^° -H^J='^ -^2;^'- (^-^"^ 

We will need the parameters tj later in the computation of the Fourier coefficients of the p function and the 
perturbation calculations for the eigenfunctions. From equation (j4.9p we know that they arc related to the 
ej as follows 



for j — 1,2. Therefore if we pick e^, we pick tj and vice versa. 

Since any potential q{x, e, a) is always Floquet isospectral to q{±x + a, e, a) where a is arbitrary, we 
cannot hope to remove the sign or translation degeneracy. We know that when e = eg that di ■ a = ui and 
S2 ■ a = 1/2, so for simplicity wc fix a so when e ~ eo then a — 0. As a result we have that 

qj{s,a,eo) = pj(s + ^,7"j) = ^ a^ cos(27rns) 

nGN 

for j = 1, 2, where the coeflficients a^ are given by Appendix \K\ We consider our manifold M{e) of potentials 
which have translation fixed as above. 

In order to prove that Af(e) actually is an analytic manifold with coordinates a = {ctj^miq)} we must 
first remind the reader of a few definitions involved in the selection of the coordinates {aj^m} defined by 
(j2.5|) as they are related to the Dirichlet spectra Hj.m{q) of the operator. We define the discriminant A(A) 
as follows 

A2(A) - 4 = 4(Ao - A) n ^^i-^l^f-^\ (4.12) 

n— 1 

Let ^rn{s,qj) = ^j^m{e,s,a) be the the solution to the system (where here we are suppressing the j) 



dfi^{e,s,a) ^^2^2 v/A2(/i„,)-4 

ds n (/^n(e,s, a) -/im(e, s,a))/n27r2 

nel, 

with ^j^m(e,0, a) — firniO,qj),k G /, where the choice of signs is initially by the sign of numerator, and 
changes whenever fXj^mi^, s, a) hits A^„j. The proof of analyticity of /i by examining (|4.13|) remains almost 
exactly the same as in [3] and is omitted here. Since there are a finite number of coordinates, it is easy to 
see that analyticity in each coordinate is preserved, and hence A/(e) is still an analytic manifold 

By McKean-Van Moerbeke [7], the initial value the sum of the initial values, /J,j,m(e, 0,a), is related to 
each directional potential qj{e,0,a) in the following way 

9j(e, 0, a) = Ao + ^ (A+„, + A^",„ - 2^j,™(e, 0, a)) 

■melj 

and this relationship remains true when the parameter s is varied 

q,ie, s, a) = Ao + ^ (A+„ + Xj^^ - 2/i,- „(e, s, a)). (4.14) 

rn (E Ij 

Using a combination of formulas on pp. 325 and 329, in [9], the eigenfunctions for each directional potential 
corresponding to XJ,^ for all j can be written as 

where Am is the zero of 4A lying between A^ and A^. It is important to note here that the formula in [3] 
is a misprint. We will also need the derivatives of the eigenfunctions which from equation (j4.15p are 

20+(e,.,a)^^-(^'^'") = y -\ (^!h(5lll^) TT A+-M.(e,.,a) 

with the derivative for (j)~ (eo, s,a) computed similarly. Let us start by considering the eigenfunctions for 
those directional potentials with j > 3. Because we are looking for the root between A^„j and XJ^^ when 
e = eo, we make the substitution A ~ ^Jm + ^ilj.m^ into (|4.12p to find that 

A2(A)-4 = e2A(l-A)/(e4A,e4) 



m(e,s,a)) = 11 -f r (4.15) 



where f{z,e4) is analytic and /(O, 0) = 7,^„ 7^ 0. Therefore for 64 sufficiently small, A„i corresponds to the 
root of 

= (1 - 2A)/(e4A, £4) + e4A(l - A)|^(e4A, £4) 

oz 



near A = i. As a result, the following estimate holds 

A^(e) - Am(e) 
giving that 



0(e4). (4.17) 



X^ie^o) ^.rnieo as) ^ 2cos^(5„(.,a)). (4.18) 

Am(eo) ^ Am(eo) 
The variable amis,a) denotes the solution to the system (|4.13p where e = eg with initial condition a under 
the change of variables (j2.5p . The same estimates above are true for the eigenfunctions 03 „j(eo7 s, a), m in 
I3 when expanded with respect to £3. We can conclude for all j > 3 

'?^iim(^o,s,Q;) = \/2cosaj^rnis,a) (4.19) 

where we know we have picked the right sign by verifying the derivative (|4.16p in the limit. 

Now we consider the case when j < 2. When e = eoj we have for all n > 1 that A^„ = A^ = ^„ = A„ so 
that terms in the product (|4.15p where n ^ m and n > 1 become 



•^jim(eo) - Aj-,n(eo) 

and for n = 1 we have 



1, (4.20) 



^i,m(<^o) - Aj,i(eo) Aj-„ - Xjs{ej) 



(4.21) 



Combining equations (|4.18p (which is still true for j < 2) and (|4.20p . we see that for e = eg, and (jjm) in 

((/.+ „(eo,a,s))^ = 2cos^(a™(s,a)) -^^ '. r — . (4.22) 

V A+„(e,)-A,-i(e,) ) 

Comparing with the derivative computed in (|4.16p we know that the correct choice of sign is 

('/'j^m(<^o,a, s)) = y2cos(a,„(s,a)) 



^tm(^j) ~ ^,,i('^j) ~ ej7i,isin^(ai,i(S'a)) 



(4.23) 



\ X\^{tj) - \j.x{tj) 

The introduction of this setup provides the necessary background to introduce the following theorem: 

Theorem 4. For all but an analytic set 0/(63, £4) in [0, 1]^, there is an open set of potentials satisfying the 
hypotheses (1),(2) and (*) in M{e) which are isospectral to only a finite number of other analytic potentials. 



In order to find the Jacobian corresponding to the invariants as given by equation (|4.7p . we must first 
figure out what it means to calculate their derivatives with respect to {aj^m} with {j,m) in Eq. We start 
with the following lemma 

Lemma 1. For {j,m) in Eq, we have 

%i(fl^^l, and 9a, is,a) ^^ when (r, fc) ^ (j, fc) 



Proof. Examining (j4.13p under the change of variables given by (|2.5p for (j, 771) in Ei and e = eo 



da 



J.™ 



(s, a) yi^],rn ~ ^o)(Aj „j - Aj- i)(Aj- ,„ - Aj- ;^) 



ds A+„ - A^ 1 - ej7ju sin ajsis, a) 



(4.24) 



Therefore aj^mis,a) depends only on aj^i and the initial data for aj^m(0,a) = Qij^m so the result follows. 
The case whenever j > 3 and e = eo, is much easier to compute. We have for all such corresponding m 

— ■'' = miT (4.25) 

ds 

so again the result follows by the same reasoning above. D 

For the computations done in the appendix, we need to know that when e^ ~ 0, (|4.23p agrees with the 
limit one would expect. In other words for (j, m) in Ei, we have 

4>tm{<^o, a, s) = ^/2cos(7^ms + a^- „,) + 0{ej) (4.26) 

which is easily verifiable by Lemma[T] and the estimates (|4.17p and (|4.20p . We have computed the cigcnfunc- 
tions in (|4.23p to illustrate that they are expressed in terms of elliptic functions, and therefore the invariants 
will not be explicitly computable. 

We can now prove the main Lemma. If we consider a potential (/(e, x, a) in M{e) then it is associated to 
a fixed set of coordinates a. Let det(J)(e,a) be the Jacobian determinant of the invariants $j^m(e,a) with 
respect to the coordinates {ctj.m} with j^m in i?Q, and det(J)(e,a) is an {N — 2) x {N — 2) determinant. 

The proof of Theorem 2] will be based on the following Lemma: 

Lemma 2. There is a choice of 61,62 in [0, 1] such that on a dense open set of a, 

det(J)(e,a) 7^0 (4.27) 

Proof. We will proceed by showing that for all fc = 1 to n — 1 

9'=det(J) 



94 



-(eo,a) = 



while 



a"det(J) 

-(60, a) 7^0 



de^ 



where n = |/i| + I/2I — 2 = |i?i|. The desired result will follows since we notice that if for some n 

d" det( J) 



de^ 



-(eo,a)7^0 and det(J)(e,a) = 



then this is a contradiction since all of the derivatives of det( J)(e, a) evaluated at any e should be identically 
zero as well, since det( J)(e, a) is an analytic function of e. 

Now we proceed to calculate the derivatives of det(J)(e,a). Let the columns Vi(6,a) of det(J)(e, a) be 
indexed by i where i ranges from 1 to iV — 2. Each i corresponds to a pair of indices (j, m) such that 

Wi(e,a) = Va$-,-^„i(e, a) 

where we are considering the pairs (j, m) ordered first by the j and then by the m. The perturbation 
calculations to find the derivatives of the invariants are located in Appendices. In order to examine the 
Jacobian further, we need the following key observations: 



dq. 



1. ^{eo,Sj • a;, a) = V(/, k) e E^^, and Vj 



2. 



dai 



-(eg, ^j • a;, a) = V(/, fc), (j, m) e Eq unless (j, fc) = (/, m) 



3. ^(eo, ^j ■ x, a) = 3," (eo, <5j • x, a) = Vj ^ 3 



The first two observations follow from Lemma [1] and formulae ()4.14p and (j4.15p . respectively. The last 
observation follows from the parametrization of the open gaps since only q^ie, S3 ■ x, a) and (/)3.m(e, S3 ■ x, a) 
for m in I3 depend on 63. 

Going back to equation ()4.7p . each invariant has the form as follows 



$j,m(e, a) 



E 

leN 
1^3 



Si ■ dj 



qi{e,Si ■ x,a) 



tmi<^^^J ■x,a)fdx 



i^j.my^ J 



(4.28) 



Now we let D denote a generic constant independent of the coordinates. When e = eo the form of the 
invariants ()4.28p for j > 3 coincides with that of [3]. Since Si and S2 form a basis for §, we know that there 
exists a nonzero pair of integers (pi^ri) such that for any third vector Si ^ S\^S-2, we have Si = piSi + riS2- 
Therefore when j > 3 



1 1 



«T2 



ITi 



*j\m(eo,a) = £» / ip2it+^-,T2))ipi{s + —-,Ti))coa^{TTm{pjS + r.it) + aj^rn)dsdt + D (4.29) 







Exactly as in [3] , we have that when ( j, m) is such that j > 3 

*&i,m(eo, a) = ci^2,jalnp^al^rj cos2aj^,n + D 

The coefficients ci^2jamp "^mr ^'^^ independent of the coordinates and nonzero. They can be found in 
Appendix VK\ However for j in {1,2}, we come across the degeneracy that 



a$, 



9a, 



l.k 



-(eo,a) = 



(4.30) 



for all {l,k) in E^. Wc know from our observations (1) and (2) that (|4.30p holds except for possibly when 
{l,k) = (j, m). In this case since again Si and S2 form a basis for S we can write 



a$ 



J,™ 



da 



(eo,a) = 



J,™ 
1 



Si-d 



-qi{eo,Si ■ x,a) 



"VTi 771/ / r \\ J 

-— (eo, Oj • X, a)) dx 



(4.31) 



-D I pi{s+—,Ti)ds —^ {eo,t,a)dt 



where I 7^ j and / is in {1,2}. But since we consider our eigenfunctions as normalized for all (j,m), e.g. 
W't'tmi^^^j ■ ^''^)II-L2(R) — I7 the right hand side of (|4.3ip is just zero. 
Therefore for all i from 1 to n we have 

Vi{eo,a) = 0. 
while for all i from n + 1 to {N — 2) we see that 



l = l,..,i-l 

(u,;(eo,a))* = ^ cx^2,i,a\p^al^r, sin2aj,™ / = i 

l>i 



(4.32) 
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Because the determinant is a multi-linear function of its rows, we may write 

dct(J)(eo,a) = det (wi, U2, ..., w„, w„+i, .., wjv-2) 
It is now clear that for all A; = 1 to 7^ — 1 

a'=det(J) 



d4 



-(eo,a) = 



however for k = n wc have 



9"det(J) ^/ NJi/^^^i ^^2 dvn 

— (eo,a) = C(n)det -— , t;— ,•■, -^—,Vn+i,...,VN-2 



del 

where C{n) is a constant depending on n only. 
From observations (1-3) we know for j in {1, 2} 



de^daij 



des '^63' ' de3 



-(eo,a)=0 



(4.33) 



except for possibly when Z = 3 or (Z,fc) = {j,m). We then note that corresponding rows with 1 < z < n in 
(|434)) take the form 



dvi 







/ = l,..,i-l 



a^*, 



= < 



3/i 



dc.,,^de, i^O,a) 1 = 1 

r > I > i 

t; — ^73— (en, a) i = r...k 

/ > r 



(4.34) 



Here the index r corresponds to (3,1) and k — r = I/3I. We can conclude from (|4.32p and (|4.34p the 
determinant (|4.33p is an upper triangular one. The determinant (|4.33p looks like 

A B 
C 

where ^ is an ?i x n block diagonal matrix, and C is an {N — 71 — 2) x {N — n — 2) block diagonal matrix. 
If the diagonal entries in the upper triangular determinant ()4.34p are nonzero, then we will arrive at the 
desired result that 



c)" det( J) 



(eo,a)^0 



(4.35) 



The collection of diagonal entries for (j, m) in Ei corresponding the block A^ for 1 < i < n are ^ — a^'"* (co; Q?). 
From Appendix [B] we know that there is a choice of ei and 62 so that these invariants are nonzero except 
on an analytic set of otj,m- Also from Appendix IB] and equation (jB.22p . whenever i > n we have diagonal 
entries corresponding to {j,m) with j > 3, corresponding to the block C are 



da-i, 



-(eo,a) 



-2c 



l,2j 



'^mpj '^mrj ^^^ ^Otj,m 



(4.36) 



These entries are only zero whenever aj,m = mod tt/2 for j > 3. The lemma is finished. 

Remark : It should be possible to remove the assumption (*) by using the standard perturbation series 
to calculate ((/)|^„ (cj , s, a))'^ around ej = 0. If S3 were generically of the from psSi + ^3(52, then wc conjecture 
that (jB.SP is nonzero provided wc expanded the eigenfunctions to order n with n satisfying the relation 
m±l = np3 or m±l = nr^ for some I in N. The calculations required to do so are difficult. This conjecture 
is discussed further in Appendix |B] D 



11 



Proof of TheoremVR This proof is very similar to the one in [3] and is again included for completeness. Let 
us start by assuming the matrix J is invertible on M(e) except for on an analytic set, say U , of (ea, 64) Recall 
that on the manifold ej and the corresponding aj^i for j = 1,2 are fixed. Then given some e with variable 
components (63,64) in [0, l]'^/U, wc let 

Since 

$(e,a) :M(e) ^ M^^^ 

the corollary follows if wc can show that the set $^^($(_F)^) is open and dense. We know the set is open 
since $~^ is open, and F is compact. If we assume that it is not dense, then the set contains contains 
an open set O which also contains a point ag which is not in F. Because the Jacobian is nonzero, $ is a 
homeomorphism on a neighborhood of ao, which implies ^{F) contains an open set. The last statement 
contradicts Sard's theorem. Now we assume that $(ai) is not in $(F) and $~^($(F)) is infinite. Let a2 
be an accumulation point of $~^(<f>(ai)). Because $ is continuous, $(012) = $(ai) and ^^ 7^ 0. It follows 
that there is a neighborhood, A'^, of a2 such that a is in N and ^(a) = $(02) implies a = a2- This is a 
contradiction to our assumption so we know $^^($(ai)) is finite. Because $ is a spectral invariant, then 
$~^($(i^)^) is a subset of the manifold which satisfies the conditions of Theorem|4l D 

This theorem has a nice corollary if wc make the following observations: 

1. Any two directions 61 and S2 form a basis for the lattice L, so our choice of basis and translate is 
arbitrary. 

2. The potentials on M{e) satisfying the conditions of the theorem are dense in the set of all analytic 
potentials in the C°° topology. 

3. The set of smooth periodic potentials which are a sum of only a finite number of directional potentials 
each with a finite number of gaps in each direction are dense in the set of finite gap periodic potentials 
in the C°°{R^/L) topology. 

4. The set of finite gap potentials is dense in the set of all C®(]R^/L) potentials in the C°° topology. 

Corollary 1. The set of analytically rigid potentials is dense in the set of smooth potentials on M^/L in the 
C°°(]RVL) topology 

A Fourier Coefficients of the p Function 

As detailed in section 1, the p- function depends on a parameter tj > 0. The complex valued function p{z, r) 
is given by 

1 1 



P(z,t) = — + J2 



z^ ^-^ \(z — n — imr)^ (n + imT)^ 

which as before is real on the line x -\ — j- and setting, 

a = e-^'^^' b = e2"(^+^) 
gives 



11 Y^ ab o\^ '^"" 

^ ^ n— — 00 ^ ' n— 1 
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Because 



and 



a" 6 



(1 - a™5)2 



Y^ w(a'"&)" "^ > 



(1 - a™6)^ 



= ^n(a-™6-i)" m<0 

n=l 



the representation 
1 



(27ri)2 
12+^ 



pix,T) 



OO CXD 



na^ e 



27rinx 



V V n(a"(™+4)e^""'' + a"('""3)e"^"'') - 2 V 



n— 1 771—1 n— 1 

Changing the order of summation we get 

1 1 °° r, — °° n 

— 1 1 ■r-^ 2na2 „ v^ "i^o. 

-r^p{x,T) = —- + y cos(27rna;) — 2 > . 

47r2 * ^ ' ^ 12 ^^ 1 - a" ^ ^ ^-^ 1 - a" 

n— 1 n— 1 

Therefore the Fourier coefficients for the p functions in the first three directions are given by: 

J _ - 871-2 nexp(-7mTj) 

" 1 — cxp(— 27rnrj) ~ 

— ncxp(— TrnTjj 



ao 






exp(— TrnTj) 



where j = 1,2. The appropriate tj will depend on the choice of ej as given in section 1. 

B Calculation of the Invariants 

In order to prove Lemma [2] we need to show that there exist ei and £2 hi [0, 1] such that 






-(eo,a) and 






(eo,a) 



are nonzero except perhaps on an analytic set of a. 
We know by (|I7)) 



$j;m(e, a) 



E A .^. gfe('^''^fc '^'O^) 



j#fc 



(5fc • dj 



f-j,m 



(e, (5j • X, a))^ (ia: 



(A.l) 
(A.2) 



(B.l) 



(B.2) 



Each qj{e, Sj ■ x, a) is independent of ea when j ^ 3. Furthermore since qk{e, Sk ■ x, a) and (^^ „(e, (5^ • cc, a))^ 
arc independent of /^j^m(e, ^j • x, a) for all j 7^ fc, so only the function (0^,j^(e, (5j • x, a))^ depends on a^^m in 
the above integral. As a result we can write 






_d_ 
des 






4 • dj 



d 

75 ('^tm(e''5j ■x,a)fdx 



(B.3) 



13 



Whenever e = eq, then g3(eo:'^3 ■ x,a) = and the derivative de^^qsieo^S^ ■ x,a) can be calculated using 
the Frcdholm alternative as in [3]- Following Appendix I of [3], we may write 

^(eo, S^- x,a) ^ V" 73. „ cos(27r53 • x + 2a3,„). (B.4) 

063 ^ 

Also from the derivation of equation (|4.19p . we can conclude that 

i(f>t,mi'^o,s,a))^ = 2cos2(7rm(,5j • x) + a,^^) + 0{ej) (B.5) 

where by Lenima[T]the order terms are bounded by ejC where C depends only on aj.,„. Hence from analytic 
perturbation theory and the derivation of (|4.19p we can use (jB.5|) to conclude that 

d(<f>+ )2 
' ^'"^ (eo.'^j • x,a) = -2sin(27r(5j • x)m + 2aj,m) + Ofe) (B.6) 

where the 0{ej) terms are bounded by e^C with C a constant depending only on the coordinate ctj^m- 
Because any two directions Si and S2 in S form a basis, we know that there exists a nonzero pair of integers 
(p;, ri) such that for any third vector Si ^ Si, Si we have Si = piSi + riSi- For easier computations we make 
the initial variable change Si ■ x ~ s and S2 ■ x = t, with the associated Jacobian, Vol(r), and rewrite the 
invariants. We also let D denote a generic constant which is independent of the coordinates, and we let 

From statements (1-3) in Section 3, (|B.4p . (|B.5|) and (|B.3p . when e — eg, we have 

(c3.,,Vol(r))-i/|i^(eo,a) = (B.8) 



4 / / I ^73,nCOs(2^n(s + t) + 2a3,n)] Pift + »j,^i)) ^;l''" {to , s , a)f ds dt 
^"e/3 / 

2 ^ 73,«aL / cos(27rn5 + 2a3,„) -^^"^ {eQ,s,a)ds 

where < j, / < 2, j 7^ /. 

When j = 1, by the hypothesis (*) on the number of open gaps that q^ has, the right hand side of (JB 
is just 

2a2^73,„ sin(2a3,™ - 2ai,„0 + 0(ei) (B.9) 

Here the C(ei) terms are bounded by eiC where the constant depends only on ai^m and a^^ for all n G I3. 
We recall that ajj — >■ as e; — >■ for all n in N and / = 1, 2 since ajj is related to ei by Equation (|4.11l) and 
(|A.ip However, we can make the constant uniform in £2- If we let 



(eo,s,a) 



Mm. < CX) (B.IO) 



sup 
s«E[o,i] Cai,ji 

then this follows from the rough estimate 
^ / / y^ 73,n cos(27rn(s + t) + 2a3^n) \ Pijt + i-^,T2)) „ ^'"^ (ep, s, a))^ - sin(27rms + 2ai.m) dsdi 

nG-fa V"6^3 / \ '"* / 

(B.ll) 

Y, ^3,nal cos(2a3.„) (Af™ + 2) < 2n (A/„ + 2) 
ne/3 
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< 



since the gap lengths ^z,n and the Fourier coefficients a^ are exponentially decreasing. Now let /3 in (0, 1) 
be a small fixed parameter. We consider the set of a such that 



|2a3^m — 2ai^m ^ kii] > /? Vfc G Z, m G /i 



(B.12) 



We let this set be denoted as Ai, and note that its complement is an analytic set. Therefore provided we 
chose ei and £2 which satisfy the inequality 



{Mm + 2)ei < 



|Qml73,^ 

2n 



■ sin(;3) 



(B.13) 



for all ra in /i and a \n Ai then (jB.Sp is nonzero for j ~ 1 and all m \n Ii. The tricky step is to prove that 
we can pick ei, £2 in (0, 1) such that IB. 131 holds for all m in /i but also so 



a2$2 



de3da2,n 



:eo,«)^0 



(B.14) 



for all n in I2 except on an analytic set of a. 

Because for small ei, a^^ > a^^ whenever ^2 > ni the right hand side of ()B.8P is already written in 
ascending order in ei for j = 2, ^ = 1. Let 



+ ^2 



(eo,a) = / cos(27rns + 2a3,„) 



da-j 



■{eo,s,a)ds. 



(B.15) 



Since we do not know if &2,m,n(eoja) = in a for all m y^ n, we pick ei as follows. Say 62,m,i(eoi Q^) is 
nonzero except on an analytic set of a, and then let the set where &2,m4(eo, a) = be denoted as ^2,m.i- If 
we can prove that for j ~ 2,1 = 1, (|B.8|) is nonzero for some a, then it will be nonzero on some open dense 
set of a's. The easiest a to select is the one when 62,m,i(eoi Q^) is at its maximum. Hence we then pick ei 
such that 



max |73,iai[62,m,i(eo,a)| > 

"6^2. m 1 



^ 73,fcai^2,rn,fc(eo,a) 



fce/3 



(B.16) 



where the max is taken over the possible values of 62,m,i(eoj <^) with a in A2_m,i, and we consider the right 
hand side of (JB.16P to be evaluated at this a as well. If b2,m..ii^0iCt) = in a, but 62. m, 2(^07 a) is nonzero 
except on an analytic set of a2,Tn, and let the set where &2.m,2(eo7<^) = be denoted as ^3 „j 2 then pick ei 
such that 



max I73, 202^2. m,2(eo,a)| > 

C«6A2,m,2 



^ 73,/cafc&2,m,A;(eo,a) 



fce/3 

fc>2 



(B.17) 



where again the max is taken over the possible values of &2,m,2(eoi Q^) with a in A2^m,2- We continue this 
process inductively. As before, let /? be a small parameter in (0, 1). We now also consider the set of a such 
that 



|2a3^„i — 2a2,jn — kn] > (3 Vfc e I^, m e h 
and let this set be denoted by A2^m,m- We know 

&2,m,m(eo,a) = sin(2a3,„i - 2q2,„i) + 0(62) 



(B.18) 



(B.19) 
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where the 0(e2) terms are bounded by e2C where is a constant depending only on a2,m and az,n for all n 
in I^. Hence our selection process terminates because 62,m.m(£07a) is not zero for a in A2^m.m provided we 
chose £2 such that 



e2\C\ < sin(/?) 



Hence wc pick ei in terms of £2 so that 



mm max 

n a£A2,m. 



(l73,na^^2,m,n(eo,a)|) 



> 



^I3,kalb2,m.k{'^0,0!.) 






(B.20) 



(B.21) 



for all m in I2 where the min„ is taken over those indices n for which 62.m,n(eoi a) is not identically zero 
in a. This choice of ei and £2 is not in contradiction to our choice of £1 small compared to £2 since the 
right hand side of the inequality (|B.21|) always has a higher order function of £1 than the left hand side. 
Furthermore 62,m,n = for all m ^ n whenever £2 = 0, so the right hand side is bounded. We conjecture 
using a computer and the standard perturbation series for 6j,rn,n(eoj ck) that the assumption qi, (72 and q^ 
have the same number of gaps could be removed. However, this is computationally difficult since it has been 
verified 6j_m,n(eo, ct) is 0{ej ) for all m up to some sufficiently large values of m and n. 

For the case with j > 3, the invariants are computed almost exactly the same way as in [3] because the 
form of the invariants coincides for these indices. In this case we have that 



*j,m(eo,a) = ci,2ja^ 



3 "'mpj ""mrj 



a„ 



cos{2aj 



D 



(B.22) 
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